Spectroscopic Evidence for Multiple Order Parameter Components in the Heavy 

Fermion Superconductor CeCoIn 5 
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Point-contact spectroscopy was performed on single crystals of the heavy-fermion superconductor 
CeCoIns between 150 mK to 2.5 K. A pulsed measurement technique ensured minimal Joule heat- 
ing over a wide voltage range. The spectra show Andreev-reflection characteristics with multiple 
structures which depend on junction impedance. Spectral analysis using the generalized Blonder- 
Tinkham-Klapwijk formalism for rf-wave pairing revealed two coexisting order parameter compo- 
nents with amplitudes Ai = 0.95 ± 0.15 meV and A2 = 2.4 ± 0.3 meV, which evolve differently 
with temperature. These observations indicate a highly unconventional pairing mechanism, possibly 
involving multiple bands. 

PACS numbers: 74.70.Tx, 74.50+r, 74.20.Rp, 74.45. +c 



The discovery of the heavy-fermion superconductor 
CeCoIn5 has attracted widespread interest in the field of 
superconductivity . Besides having the highest critical 
temperature T c = 2.3 K among heavy-fermion materi- 
als, CeCoIns also shares some unconventional properties 
with the high-T c cuprates. First, CeCoIn 5 has shown pro- 
nounced non-Fermi liquid behaviors, suggestive of quan- 
tum critical phenomena that could arise from competing 
orders 0, 0|- Second, CeCoIns has shown low-energy 
quasiparticle excitations and a power-law temperature 
dependence in the NMR spin relaxation, indicative of 
nodes in the superconducting energy gap 0, 01 El 0- 
These nodal charateristics are consistent with e?-wave 
pairing symmetry |8j, which could be produced by an- 
tiferromagnetic fluctuations 0. Unlike the cuprates, on 
the other hand, CeCoIns is an intermetallic compound 
with multiple sheets on the Fermi surface [lfj, [llj . Such 
complex Fermi topology could involve several bands in 
the pairing pro cess, giving rise to multiple pair poten- 
tials 0EIII3. 

Point-contact spectroscopy (PCS) has been a proven 
microscopic technique for studying unconventional su- 
perconductors. For the high-T c cuprates, PCS pro- 
vided the earliest measurements of the superconducting 
gap spectra 0. In MgB2, PCS was key in revealing 
two coexisting s-wave gaps [lfij . PCS has been previ- 
ously performed on several heavy-fermion superconduc- 
tors 0, 0, 0, |2(J, E| • For superconductors with gap 
nodes, PCS can i n ge ner al p rovide information on the 
pairing symmetry |22l I23L 124. l25| . In this Letter, we re- 
port PCS measurements on single crystals of CeCoIns in 
the temperature range 150 mK to 2.5 K. We observed 
Andreev-reflection characteristics with multiple struc- 
tures, whose dependence on junction impedance indicates 
two coexisting order parameter components (OP) with 
nodal characteristics. These OP's show sizable ampli- 
tudes relative to T c and different evolutions with temper- 
ature. Our observations suggest a highly unconventional 



pairing mechanism, possibly involving multiple bands. 

In PCS, electronic transmission between a normal 
metal and a superconductor is measured as conductance 
dl/dV versus bias voltage V across a ballistic contact 
junction. For a transparent contact, dl/dV is primarily 
determined by Andreev reflection, based on the conver- 
sion of electrons or holes into Cooper pairs, which dou- 
bles dl/dV inside the superconducting energy gap. For 
non-transparent junctions, dl/dV involves both Andreev 
reflection and quasiparticle tunneling. The standard 
model for calculating dl/dV was given by the Blonder- 
Tinkham-Klapwijk (BTK) theory for s-wave pairing [261, 
and subsequently generalized for <i-wave pairing |27l |28| . 
A key spectral signature of the d-wave gap nodes is the 
zero-bias conductance peak, which arises from surface 
states bound by phase interference between consecutively 
Andreev-reflected quasiparticles j2^- This peak struc- 
ture is to be distinguished from the hump structure as- 
sociated with conventional Andreev bulk states. In the 
generalized BTK scenario, relative manifestation of the 
Andreev surface versus bulk states depends on both junc- 
tion orientation and a dimcnsionless parameter Z repre- 
senting junction impedance [30], thus allowing the OP to 
be studied jHHI. 

The single crystals of CeCoIns used in this work were 
grown by a self- flux method l], and characterized by 
both x-ray diffraction and magnetic susceptibility to con- 
firm material uniformity. The crystals were platelets ap- 
proximately 1x1x0.2 mm 3 in size, each showing a sharp 
superconducting transition at T c = 2.3 K. The crystal 
surfaces were etched with HC1 and rinsed with ethanol 
prior to measurement, in order to remove any residual In 
flux. High purity Pt-Ir tips were used as normal-metal 
electrodes, gently pressed onto the c-axis face of each 
crystal with a spring-cushioned differential micrometer. 
This point-contact mechanism was attached to the mix- 
ing chamber of a high cooling-power 3 He/ 4 He dilution 
refrigerator, and enabled the junction impedance to be 
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FIG. 1: Normalized dl /dV vs V data for Pt-Ir point contacts 
on CeCoIns at 0.43 K. Top panel (a) is for a 0.4 f2 junction. 
Bottom panel (b) is for a 0.2 Q junction, with the 1.5 K curve 
(dashed) also plotted to clearly show the double humps. 



varied in situ at low temperatures. The point contacts 
we measured were in the 0.2-1 Q range, consistent with 
the contact size being in the ballistic regime [3l| . To min- 
imize Joule heating in the junction over a wide voltage 
range, our spectroscopy data was acquired by a pulsed 
technique: 2 ms current pulses were applied through the 
contact in 20% duty cycles, and the junction voltage was 
measured 80 times within each pulse and then averaged. 
The I vs V curves were obtained by varying the current 
level, and then numerically differentiated to obtain the 
dl/dV vs V spectra. 

Two types of spectra were observed in our mea- 
surements, depending on the point-contact impedance. 
These measurements were reproducible on multiple spots 
over different samples and repeated at each spot to rule 
out any surface destruction by the point contact. Fig- 
ure ^ shows dl/dV spectra taken at 0.43 K well below 
T c , after normalization relative to spectra taken above 
T c . The top panel is for a 0.4 SI junction, and the bottom 
panel is for a 0.2 SI junction. Distinct spectral features 
are seen in the top panel, with a sharp zero-bias peak 
dipping at ~ ± 1 mV into a broad spectral hump ~ ± 
2.5 mV in width. Small kinks are also visible on the peak 
at ~ ± 0.3 mV and ~ ± 0.5 mV. The main peak, dip 
and hump structures evolve differently with decreasing 
junction impedance. As seen in the bottom panel, the 
peak becomes an asymmetric inner hump ~ ± 1 mV in 
width, the dips get filled in, while the outer hump re- 
mains largely unchanged. These hump structures are the 
classic signatures of Andreev reflection, which introduces 
excess spectral states inside the energy gap |2(|. These 
excess states expectedly diminish with temperature, as 



is evident in the 1.5K data [dashed curve in Fig. 1(b)]. 
The zero-bias peak, on the other hand, is key evidence 
for nodes in the gap j^. It is worth noting that peak 
and hump structures of similar shapes and energy scales 
have been reported in an earlier PCS study of CeCoIns , 
although appearing separately in different spectra |32| . 
Our measured spectra are clearly hybrid in character, 
each containing multiple structures. 

To identify the multiple spectral features observed in 
our data, we consider theoretical spectra from the gen- 
eralized BTK model. Shown in the top panels of Fig. El 
are the simulated dl/dV spectra for a d-wave OP, plot- 
ted in normalized units vs eV/A, where A is the d- 
wave gap maximum |33j | . The choice of ci-wave symmetry 
here is motivated by both thermodynamic and transport 
data 0, , and intended to illustrate the generic spectral 
dependence on junction orientation and impedance. The 
curves in Fig. 2(a) are for a high-impedance (Z = 1) 
junction, and the curves in Fig. 2(b) are for a low- 
impedance (Z = 0.5) junction. In each plot, the dot- 
ted/solid curve is for a nodal/antinodal junction (nor- 
mal to a nodal/antinodal axis), while the dashed curve 
models the effect of junction roughness by averaging over 
all intermediate orientations. Note that an ideal c-axis 
junction would produce similar spectra as the antinodal 
case, since there is no OP sign change about the junction 
normal in either case to allow for Andreev interference. 
The overall spectral evolution between peak and hump 
structures is a direct manifestation of the competition 
between Andreev surface and bulk states |27j . 

From these generic spectral simulations, the data in 
Fig. H can be interpreted as the superposition of two 
types of spectral contributions. Namely, the sharp peak 
structure comes from Andreev surface states due to a 
high-Z nodal junction, and the broad hump structures 
come from Andreev bulk states due to low-Z antinodal 
junctions. The appearance of two effective Z's, with 
very different dependences on junction impedance, is in- 
dicative of different Andreev coupling to two distinct 
OP's. To demonstrate this two-OP scenario, we have 
developed a superposition model, based on the "serial" 
precedence of surface over bulk states in junction trans- 
mission. More specifically, when bulk spectra from two 
different OP's coexist, their superposition is essentially 
additive 0], since bulk states can be accessed in "par- 
allel". However, when both surface and bulk spectra 
are involved, the junction transmission becomes effec- 
tively "serial", thus justifying a multiplicative superpo- 
sition within energies (\eV\ < A) where Andreev sur- 
face states can readily form. This serial model is demon- 
strated in Fig.[2fc), by superposing a peak spectrum (left 
inset) with a hump spectrum (right inset) of triple the en- 
ergy scale (i.e. A 2 = 3Ai). Here the component spectra 
were multiplied for \eV\ < Ai and added for \eV\ > A\, 
following our model justifications. The peak-dip-hump 
structures seen in the data of Fig. ^a) are remarkably 
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FIG. 2: Spectral simulations using the generalized d-wave 
BTK formalism, (a) and (b) are for Z = 1 and Z = 0.5 junc- 
tions, with nodal (dotted curves), antinodal or c-axis (solid 
curves), and angle- averaged (dashed curves) orientations, (c) 
and (d) show serial and parallel superpositions of two spec- 
tral contributions (insets) for two OP's with A2=3Ai and 
different Z's. 



well reproduced here in Fig. EJc). For comparison, the 
parallel model superposing two low-.Z bulk spectra (in- 
sets) is shown in Fig. 2(d), also generically reproducing 
the multiple- hump data seen in Fig. IHb). The overall 
spectral resemblance between our simulations and data 
is robust evidence for the coexistence of two OP's. 

Some general remarks about our two-OP spectral anal- 
ysis should be made. First, our model was intended 
to show generically how two coexisting OP's with gap 
nodes could produce the multiple spectral structures ob- 
served. The distinctively serial relationship between the 
peak and hump structures clearly establishes the surface- 
state nature of the former, as arising from Andreev inter- 
ference for a nodal OP. However, although our data can 
be explained within a d-wave framework, we cannot rule 
out the presence of other OP line or point nodes, along 
either the pole or the equator, such as in the case of 
UPt3 34]. Precise determination of the pairing symme- 
try in CcCoIn5 would require a systematic study of the 
spectral anisotropy |2^| , along with an extension of the 
generalized BTK theory beyond its two-dimensional for- 
mulation. Second, the non-trivial spectral evolution we 
observed versus junction impedance indicates a complex 
fc-space dependence of Z , with the nodal-junction states 
dominating at high Z and antinodal-junction states dom- 
inating at low Z. While surface roughness could allow 
for nodal-junction surfaces to exist on a nominally o 
axis crystal, a detailed explanation of the peak-to-hump 



evolution would require full understanding of how Z de- 
pends on the complex band structure of CeCoIn 5 pjof. 
For example, multiband coupling could in theory affect 
the formation of Andreev surface states [35J. The effects 
of band structure on quasiparticle tunneling have also 
been studied j3||- Third, the spectral heights tend to be 
smaller in the data than in the model, a difference which 
could be attributed to non-superconducting spectral con- 
tributions from either uncondensed quasiparticles j3^] or 
Kondo scattering pjjj ]. 

The temperature dependence of our spectral data was 
also examined. Figure shows spectral evolution of the 
data from Fig.^in the temperature range 150 mK to 2.5 
K. The spectra are staggered for clarity, with arrows in 
Fig. 3(a) to indicate the two-OP amplitudes determined 
from the serial model above, and a dotted baseline in 
Fig. 3(b) to underscore the "excess" spectral area asso- 
ciated with bulk Andreev states. The OP amplitudes 
Ai(T) and A 2 (T) are plotted in Fig. gfc), along with 
theoretical (dotted) curves calculated from the BCS gap 
equation. The excess spectral area S is similarly plotted 
in Fig. ISfd), after normalization by its base-temperature 
value Sq [ljj. From Fig. 3(c), it is clear that both OP 
amplitudes approach distinct zero-temperature values, 
Ax = 0.95±0.15 meV and A 2 = 2.4±0.3 meV, and vanish 
near T c — 2.3 K, consistent with both OP's being com- 
ponents of the same superconducting order. This com- 
mon T c also argues against the presence of a proximity- 
induced superconducting layer in our junctions, which 
should cause the smaller order parameter component to 
vanish below the bulk T c [lfij |. However, while Ai(T) 
is well described by the BCS gap equation, A 2 (T) de- 
viates markedly from mean-field behavior. This devia- 
tion is also evident in the reduced spectral area (S/Sq) 
plot in Fig. 3(d), indicating a predominance of the larger 
OP for parallel superposition. Similar deviations have 
been observed in other heavy-fermion superconductors, 
and attributed to the nodality of highly complex pairing 
symmetries 

HUE El Alternatively, the difference be- 
tween Ai(T) and A 2 (T) could be the signature of novel 
interplay between two different types of order j^, . 

Finally we discuss the physical implications of our re- 
sults on the pairing mechanism in CeCoIn 5 . First, assum- 
ing that each of the two energy scales identified above can 
be directly assigned to a superconducting OP, they would 
correspond to gap-to- T c ratios of 2Ai/fcsT c = 9.5 ±1.5 
and 2A 2 /fcsT c = 24 ± 3. These ratios are much larger 
than the BCS weak-coupling value of 3.5 for phonon- 
mediated pairing, and well beyond the strong-coupling 
limit even after d-wave corrections |4l|. One conceivable 
way to enhance the gap-to- T c ratio is through inter-band 
coupling, whereby carriers from different bands could in- 
teract to result in multiple pair potentials sharing a com- 
mon T c [l2^ . This multi-band scenario would be physi- 
cally plausible for CeCoIns, considering that its Fermi 
surface has four distinct sheets with different topologies 
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FIG. 3: Temperature dependence of the data from Fig. Q 
A subset of the spectral evolutions are shown in (a) and (b). 
The OP amplitudes Ai(T) and A 2 (T) determined from (a) 
are plotted in (c). The reduced spectral area S/Sq extracted 
from (b) is plotted in (d). Theoretical BCS curves (dotted) 
are included to indicate deviations from mean-field behavior. 



and effective masses Furthermore, Andreev scat- 

tering for a heavy-mass 2D sheet would be inherently 
weaker than for a light-mass 3D sheet, due to poorer 
Fermi- velocity matching across the junction j3fj| . This 
multi-band effect could provide a natural explanation 
for the two different Z scales observed in our spectra. 
However, even allowing for inter-band coupling between 
highly disparate densities of states ^3 1 a sizable "intrin- 
sic" 2A/fcsT' c , intermediate between rs 9.5 and 24, may 
still be needed to explain our data 0, \^ . Such an in- 
trinsically large gap-to-T c ratio would present a serious 
challenge to current theoretical formulations fill |43| , at 
least within the Fermi-liquid framework, thus indicating 
a highly unconventional pairing mechanism in CeCoIns. 

In summary, we have performed point-contact spec- 
troscopy on the heavy- fermion superconductor CeCoIns. 
Andreev-reflection characteristics with multiple struc- 
tures were observed. Spectral analysis using the gener- 
alized BTK formalism revealed two coexisting order pa- 
rameter components with nodal symmetry and sizable 
amplitudes. These observations suggest a highly uncon- 
ventional pairing mechanism in a multi-band scenario. 
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